Excess Modes in the Vibrational Spectrum of Disordered Systems and the Boson Peak 
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We study a disordered vibrational model system, where the spring constants k are chosen from a distribution 
P(k) oc 1/k above a cut-off value k m i n > 0. We can motivate this distribution by the presence of free volume 
in glassy materials. We show that the model system reproduces several important features of the boson peak in 
real glasses: (i) a low-frequency excess contribution to the Debye density of states, (ii) the hump of the specific 
heat cv (T) including the power-law relation between height and position of the hump, and (iii) the transition to 
localized modes well above the boson peak frequency. 

PACS numbers: 63.50.+X, 65.60.+a 



The vibrational spectrum and the corresponding thermal 
properties of a large number of disordered materials exhibit 
two characteristic anomalies [jj],^]. While the anomalous lin- 
ear low-temperature specific heat can be described well by 
the model of two-level systems the origin of the second 
important anomaly, the boson peak, is still under discussion. 
It is not clear, if both anomalies are related [Q]. The boson 
peak refers to an excess contribution to the usual Debye den- 
sity of states (DOS) at low frequencies [^]— P~0|] - In silicate 
and germanate glasses, for example, the boson peak frequency 
t^pcak was reported between 4.5 and 13.5 THz, about 16 — 40 
times below the upper band edge loe In silicate glasses 

[pIPIJTol] and polybutadiene and polysterene glasses JTc|], the 
peaks occur at frequencies about 4 — 9 times smaller than the 
corresponding Debye frequency wd- 

The boson peak also shows up as a small hump in the scaled 
specific heat cv(T)/T 3 in the temperature range 2 — 30 K 
©EErl H ^ ne num P s f° r different materials have a com- 
mon shape. Position and height follow a power-law depen- 
dence, 
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hump 
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with x close to 2 [ |1 1| | (see Fig. 3(b)). The hump in the specific 
heat is more pronounced for strong and intermediate than for 
fragile glasses Apart from molecular dynamics simula- 
tions on relatively small systems that do not yet give a clear 
picture of the origin of the boson peak [oL relaxational and 
vibrational models have been proposed. Many experiments 
indicate that the relaxation-type soft potential model [[l3|] is 
appropriate for fragile glasses, while vibrational models are 
more appropriate for strong and intermediate glasses [ |l4] , |T5| l 
(see also [|l^|). 

In this paper, we do not wish to enter the controversial dis- 
cussion of the microscopic origin of the boson peak. Instead, 
we concentrate on an idealized disordered vibrational model 
system that can be solved by standard numerical procedures. 
We show that the model system is able to reproduce several 
important features of the boson peak in glassy materials. 

Our starting point is the same as the one of Schirmacher 
et al. JT^], who considered a simple cubic lattice of coupled 
harmonic oscillators with random scalar spring constants be- 
tween nearest-neighbor masses. The spring constants were 



chosen from a truncated Gaussian distribution with a lower 
cut-off value fc m ; n . As a result of the disorder, a hump in the 
scaled DOS Z(ui)/oj 2 occurs at relatively high frequencies. 
When shifting /c nun towards negative values (which imposes 
some unrealistic instability to all infinite lattices) the hump is 
shifted towards frequencies about 2 — 3 times smaller than wn, 
which is still far from the frequency range where the boson 
peak occurs. Model systems with homogeneous distributions 
of positive spring constants (see below) yield maxima in the 
same frequency range as the Gaussian model system. 

In this work we also consider a disordered vibrational 
model system. The scalar spring constants k are chosen from 
a power-law distribution P(k) oc 1/k above a positive cut-off 
value fc m ; n . We can motivate P(k) by assuming a distribu- 
tion of the free volume. We find that a low-frequency peak 
occurs in the DOS at w pca k, in a frequency range compara- 
ble to the one of the boson peak of real glasses. We expect 
that in our model system external pressure will shift the peak 
towards larger frequencies, an effect which has also been ob- 
served experimentally. Investigating the localization proper- 
ties of our model system we find that the modes around the 
peak are extended. Localized modes occur above a crossover 
frequency oj c « 3o; pca k- This interesting property seems to be 
in line with experiments in real glasses. The peak at w pea k is 
also visible in cy(T) as a small hump, which shows similar 
features as the boson peak in real systems, among them the 
power-law dependence (|l[). 

To be specific, we chose the spring constants from a nor- 
malized distribution function 



P(k) 



1 



1 



ln(fc max /fc m i n ) k 



k e [k, 



mm i "-■maxj ? 



(2) 



where fc max / k m [ n , the only parameter in the distribution, con- 
trols the amount of small spring constants. The form of Eq. ^ 
can be motivated as follows: Consider a disordered system, 
where the neighbor distances between the particles are 
in the interval [a m i n ,a max ]. The limits a min and a max are 
related to the size of the particles and the maximum diam- 
eter of holes in the system, respectively, and thus character- 
ize to some extent the distribution of the free volume in the 
system. For simplicity, we assume that the distances are dis- 
tributed homogeneously between a m i n and a max . When the 



1 



distances between two neighboring masses i and j fluc- 
tuate, the spring constants kij between them also fluctuate. 
We assume that the kij decay roughly exponentially with a%j, 



exp[(a n 



j)/a*] with a characteristic de- 



cay length a*, which might be a reasonable assumption in a 
strong or intermediate glass with covalent binding [ |l5| ] and 
negligible charge separations. Combining both assumptions, 
we arrive at Eq. (Q), with the control parameter 



c/ fcn 



exp[(a„ 



OK 



(3) 



Note that the distribution (||) is equivalent to an homogeneous 
distribution of the logarithm of spring constants, P(log k) = 
const. 

We have calculated numerically the DOS of this model sys- 
tem, where in a simple cubic lattice unit masses m at nearest- 
neighbor sites i and j are connected by springs kij chosen 
from P(fc). Assuming scalar coupling constants kij, the dif- 
ferent components of the displacements decouple and we ob- 
tain the same equations of motion 



/(*) 



dt 2 



= ^2 hj [ui(t) -Uj(t)\ 



(4) 



for all spatial components of Uj. The sum runs over the 
nearest neighbor sites i of site j. The ansatz Uj(t) — 
t/j" exp(—ioj a t) leads to an homogeneous system of equations 
for the N unknown ipj, from which the N real eigenvalues 



> and the corresponding eigenvectors (ipii 
- 1, . . . , N can be determined. 
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FIG. 1.: (a) Density of states (DOS) Z(uj/uj ) and (b) rescaled 
DOS Z(lo/u!q)/ (uj/ujo) 2 versus uj/ujo for our model system (ujq = 
fcmai/m). The symbols correspond to fc m in/fcmax = 10 _J (A), 
10~ 4 (O), 1CT 3 (□), and 10~ 2 (O). The numerical results for an 
homogeneous distribution fc G [0, fc max ] are also shown (dashed 
curves with + symbols). In the inset, the frequency positions 
Wpeak/wo of the maxima of Z(uj/ujo)/(lo/ujo) 2 are shown ver- 
sus the relative maximum size of holes (a max — a m i n )/a* = 

In (fcmax / fcmin) • 

To calculate the vibrational properties and the related spe- 
cific heat cy of the model system, we have employed the 



method of Williams and Maris [ |l8| ] and the Lanczos algo- 
rithm JT^]. In F ig. 1 we show the DOS Z(uj/lu ), with 
c^o = -\/k max /m, and the rescaled DOS Z (lu / uj ) / (lo / u>q) 2 
for several values of fc m in/fcmax between 10 -2 and 10 -5 [ p0| 
and system sizes of up to 65 3 masses. For all values of 
fcmin/fcmax, both quantities show a broad maximum at fre- 
quencies a;p Ca k, which become smaller, if fc m i n /fc max is de- 
creased. Well below w pea k, we observe the conventional De- 
bye behavior, Z(uS) oc uj 2 . For comparison we also calculated 
the DOS of a vibrational model system with an homogeneous 
distribution of spring constants k G [0, fc max ] where the maxi- 
mum in the DOS simply results from a broadening of the van 
Hove singularity. The systems with power-law distributed kij 
have their maxima in the low-frequency regime, where also 
the boson peak in glasses is observed. Estimating the Debye 
frequency uj o from the plateaux in the scaled DOS (Fig. 1(b)), 
we find for the maxima in Z(ui)/uj 2 the values Wrj /^ pC ak ~ 3 
for fc min = 1CT 2 , oj D /oj peak ~ 4.5 for fc min = 1(T 3 and 
wri/wpcak ~ ^ for fe m i n = 10~ 4 and still larger values for 
smaller fc mm . In contrast, the peak for the homogeneous dis- 
tribution, shows up at frequencies cj poa k only a factor of 2 
below wd ■ 

When fc m in/ fcmax is enhanced, the maximum in Z(lS) tends 
to higher frequencies. Since in our model system fc m in/fcmax 
is related to the free volume parameter (a ma x — &mm) via 
Eq. (|3|), the peak is shifted towards higher frequencies, when 
(a ma x - a m in) is reduced. We can imagine that (a ma x - a m in) 
can be reduced by applying external pressure to the system. 
Indeed, experiments on glasses indicate that the boson peak 
is shifted to higher frequencies, when the glasses were per- 
manently desified under pressure |^]. This is in line with the 
results of our model system. 
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FIG. 2.: (a) Rescaled specific heat cv(T)/T 3 versus temperature 
T for fcmin/fcmax = 3 x (□), 10~ 4 (O), 3 x 10" 4 (A), IO" 3 

(V)> 3 x 10 -3 (O), and 10 -2 (+) and for an homogeneous distri- 
bution fc G [0, fcmax] (x). In (b) the results from (a) are rescaled 
by dividing T by T hump and c v (T)/T 3 by c v (r hum p)/T 3 ump . 
Deviations from the common shape occur for the smallest and the 
largest parameter as well as for the homogeneous distribution of 
spring constants. The experimental values for several silica glasses 
(large symbols, redrawn from Fig. 1 of Wm), are also shown. 
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Next we consider the specific heat cy(T), which is related 

to Z(w) by: 

2 



cy(T) = fcfl / duiZ(ui 



f huj y 



>/t B T 



(5) 
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Figure 2(a) shows cv(T)/T 3 versus temperature T for 
several values of fc m m/fc max as well as for the homoge- 
neous distribution of spring constants. As expected from 
the behavior of Z{u>), a maximum in cy(T)/T 3 occurs at 
Thump, which is shifted towards lower temperatures and in- 
creases in height, when k m m/k m ax decreases. For test- 
ing for a common shape of the hump, we have plotted 
[cy(T)/T 3 ]/[cy(T hump )/T h 3 ump ] in Fig. 2(b) as a function 
of T/Thump and compared the results with the experimental 
curves for several vitreous silica [O]. It is interesting that 
the results for our model system at intermediate fc m i n /£: max 
values agree reasonably well with the experimental data. No 
fitting parameter was involved. The figure shows also that 
major deviations from the common shape occur for the homo- 
geneous model system in the low-temperature range. 
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FIG. 3.: (a) The values of the maxima of c v {T)/T 3 from 
Fig. 2(a) versus the positions Th ump of the maxima. In the log- 
log plot the data fall onto a straight line described by a power-law: 
cy(T hump )/T h 3 ump ~ T-J np with x » 2.1 (dotted line), (b) Ex- 
perimental results for several glasses and amorphous solids (•) as 
well as for some crystalline materials (A), redrawn after Fig. 2 of 
[ pl| |. The data for the glasses and the amorphous solids are consis- 
tent with our theoretical results x w 2.1 (dotted line). If also crys- 
talline systems are included, one obtains a smaller exponent. Note 
that the scales of both axes are only changed by constant factors in 
(a) and (b), so the range is exactly the same. 

To see, if the specific heat also satisfies the power-law re- 
lation ([!]) we have plotted cy (Thump) /T 3 ump as a function of 
Thump in a double-logarithmic fashion. The result, shown in 
Fig. 3(a), agrees surprisingly well with the experimental data 
for a large number of glassy and amorphous solids, shown in 
Fig. 3(b). The basic power-law ([!]) is reproduced with an ex- 
ponent x that is, within the error bars, identical with the one 
of real glasses. Note that both, Figs. 3(a) and 3(b), show the 
same range for the x- and y-axes. Again, no fitting parameter 
was involved. It is interesting to note that also for crystalline 
systems the same power-law behavior has been found exper- 
imentally. An explanation of this is beyond the scope of this 
paper. 



Finally we turn to the localization behavior of the vibra- 
tional modes in our model system. By general arguments, we 
can expect extended phonons for small frequencies w < w c 
and localized modes for high frequencies w > u> c . The local- 
ization mechanism for vibrations is quite similar to the Ander- 
son localization of quantum particles in disordered solids (see 
also [ pTp2| ]) and it is interesting to see, if the modes in the 
peak regime fall into the extended or into the localized range. 

To estimate the ratio between w pca k and w c , we have em- 
ployed the method of level statistics [Q. Level statistics 
(with a fixed system size) have been applied earlier to the 
Gaussian model system indicating extended modes in the 
vicinity of w pea k. In the method, one first calculates the eigen- 
values w 2 of the vibrational equation and then determines the 
eigenvalue spacings s a = (w 2 — w^_ 1 )/A, where A is the 
mean eigenvalue spacing in the frequency range considered. 
The dependence of the mean squared eigenvalue spacing (s 2 ) 
on the system size N indicates the localization behavior. With 
increasing system size, Jo = ^(s 2 ) tends to J lgner ~ 0.637 
for extended modes, while it approaches jP° lsson — l for lo- 
calized modes. Figure 4 shows the results of this analysis for 
fcmin/fcmax = 10~ 2 and 10" 3 . From the figure we can de- 
duce w c w 1.3 wo for fcmin/fcmax = 10~ 2 and W c fti 0.7 Wo 
for fcmin/fcmax = 10~ 3 . For fc min //c max = 1(T 4 we obtain 
w c « 0.4 wo in an analogous way. A comparison with the 
position of the peak (taken from the inset of Fig. 1(b)) yields 
w c « 2.7w pcak , w c w 3.1w pcak , and w c « 3.9w pcak for 
^min/fcmax = 10~ 2 , 10~ 3 , and 10~ 4 , respectively. 
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FIG. 4.: 

fcmin/fcmax = 10" 2 and (b) fc min /fc m ax = 10~ S . J {lO,N) 

is shown versus cj/cjo for different system sizes TV. Values of 
Jo « jWignm _ g ggy indicate extended modes, while Jo ta 
jPoisson _ ^ i n di ca tes localized modes. In the insets, magnifica- 
tions of the transition regions are shown, which allow to determine 
the transition frequencies ui c from the intersection of the curves. 

Experimentally, it is not settled whether the frequency 
w pca k of the boson peak is close to the localization dereal- 
ization frequency w c or well below. Scattering experiments on 
vitreous silica have been discussed controversially. One inter- 
pretation yields localized modes at the boson peak frequency 
w poak ~ 6 — 9 THz p4|], while others indicate the presence 
of propagating modes up to w c « 20 THz well above 
the boson peak. Very recent scattering experiments on glassy 
glycerol yield extended modes up to at least w c rj 2.5 w poa k 
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The results from [25 1 and [[260 are in line with the results 
of our model system. 

In summary we have considered an idealized disordered vi- 
brational model system, consisting of a simple cubic lattice of 
coupled harmonic oscillators with random scalar spring con- 
stants between neighbors. The spring constants were chosen 
from a power-law distribution P(k) oc 1/fc above a positive 
cut-off value fc m ; n . We have motivated the model by the pres- 
ence of free volume in glassy materials. Since the model de- 
scribes (surprisingly) well several non-trivial features of the 
boson peak in glassy systems, we believe that despite its sim- 
plicity it may capture essential parts of the physics of this odd 
phenomenon - in particular as no single fit parameter was in- 
volved. 

We like to thank Prof. Michael Klinger for illuminating 
discussions on the subject. 
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